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We present results for phase structure of lattice QCD with two degenerate flavors (Nf — 2) 
of Wilson quarks at finite temperature T and small baryon chemical potential /is. Using the 
imaginary chemical potential for which the fermion determinant is positive, we perform simulations 
at points where the ratios of pseudo-scalar meson mass to the vector meson mass m^-/m p are 
between 0.943(3) and 0.899(4) as well as in the quenched limit. By analytic continuation to real 
quark chemical potential fi, we obtain the transition temperature as a function of small /is. We 
attempt to determine the nature of transition at imaginary chemical potential by histogram, MC 
history, and finite size scaling. In the infinite heavy quark limit, the transition is of first order. 
At intermediate values of quark mass m q corresponding to the ratio of m„-/m p in the range from 
0.943(3) to 0.899(4) at a^ii — 0.24, the MC simulations show absence of phase transition. 

PACS numbers: 12.38.Gc, ll.10.Wx, 11.15.Ha, 12.38.Mh 



I. INTRODUCTION 

QCD at finite temperature and density is of funda- 
mental importance, both on theoretical and phcnomeno- 
logical grounds. It describes relevant features of particle 
physics in the early universe, the neutron stars and the 
heavy ion collisions. At high density and low tempera- 
ture, some QCD-inspired models suggest a complicated 
phase structure [H, and at sufficiently high temperature 
and small density, QCD predicts a transition (In this 
paper "transition" refers to the change in dynamics, ir- 
respective of the order of the phase transition.) from low 
temperature hadronic matter to high temperature quark 
gluon plasma (QGP). Probing this transition is one of the 
main purposes of the experiments of SPS, LHC(CERN) 
and RHIC(Brookhavcn). Because QCD is strongly in- 
teracting, perturbative methods do not apply, and the 
only first principles method to investigate these transi- 
tions is by means of lattice Monte Carlo (MC) simula- 
tion. However, lattice MC simulation is based on im- 
portance sampling, which can not be directly applied to 
the nonzero baryon density case because of the complex 
fermion determinant [|J for SU(3) gauge theory. 

Enormous efforts have been made to solve this complex 
action problem. Fodor and Katz used a two-dimensional 
generalization of the Glasgow reweighting method[3( to 
study the phase diagram of lattice QCD with Kogut- 
Susskind (KS) fermions [H; Allton et o/.@] attempted to 
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improve this method by Taylor expansion of the fermionic 
determinant and observables around fi = 0. 

The imaginary chemical potential method 0,(3 nas a ls° 
been employed to circumvent the "sign problem" . D'Elia 
and Lombardo.Tj applied it to investigate the phase di- 
agram of lattice QCD with four flavors of KS fermions. 
De Forcrand and Philipscn studied the phase diagram of 
lattice QCD with two flavors [6(, three flavors and (2+1) 
flavors 8] of KS fermions. 

Monte Carlo simulation with imaginary fi has a cou- 
ple of technical advantages. It is computationally simple 
and allows control over the systematic error by fitting the 
nonperturbative data to a Taylor series. Furthermore, it 
is a good testing ground for effective QCD models: ana- 
lytic results can always be continued to imaginary \i and 
be compared with the numerics there. The main disad- 
vantage of this approach is its limitation to the range 

H/r<7r/30. 

The KS fermion and Wilson fermion approach have 
their own advantages and disadvantages. The KS fermion 
formalism preserves the U(l) chiral symmetry, whereas it 
does not completely solve the species doubling problem. 
One staggered flavor at lattice corresponds to four flavors 
in the continuum limit and in simulation the fermion de- 
terminant is replaced by its fourth root. Such a replace- 
ment is mathematically unjustified 9], and it might lead 
to the locality problem in numerical simulations [lOj. In 
Ref. [ll[ , it is pointed out that the fourth root of the stag- 
gered fermion determinant has phase ambiguities which 
become acute when Ke(fi) exceeds half of the pion mass. 

Although Wilson fermions explicitly break the chiral 
symmetry which is one of the most important symme- 
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tries of QCD, they completely solve the species doubling 
problem. So it is of interest to investigate QCD phase 
diagram with them. 

In this paper, we attempt to investigate lattice QCD 
with two degenerate flavors of Wilson fermions. In 
Sec. UH we define the lattice action with imaginary chem- 
ical potential and the physical observables we calculate. 
Our simulation results are presented in Sec. IIIII followed 
by discussions in Sec. IIVI 



The Polyakov loop L is defined as the following: 



II. LATTICE FORMULATION WITH 
IMAGINARY CHEMICAL POTENTIAL 

The partition function of the system with Nf degen- 
erate flavors of quarks with chemical potential on the 
lattice is 



(1) 



Z = I [dI7][<ty][#]e _s «- s ' 
= J[dU] (BetM[U]) Nf e- s " 



where S g is the Yang-Mills action, and Sf is the quark 
action with the quark chemical potential /i. Here /x = 
Hr + ifij, )j,r, /i/ £ 1Z. For S g , we use the standard 
one-plaquette action 



s 9 = -fE Tr ^+ t, p t - 2 )> 



(2) 



where (3 — 6/g 2 , and the plaquette variable U p is the 
ordered product of link variables U around an elementary 
plaquette. For Sf, we use the the Wilson action 

f=l x,y 

where re is the hopping parameter, related to the bare 
quark mass m and lattice spacing a by re = l/(2am + 8). 
The fermion matrix is 



M Xjy (U, re, fi) = S XtV - K,y~] 

4=1 



(1 - 74)e a "P4W^, s _4 



+ (I + 7 4)e-^C/ 4 t (x-4)^ ia+ 4 



(4) 



In this paper, we use as our observables the mean value 
of the plaquette which we denote by P, the Polyakov loop 
L and the chiral condensate tpip, we also calculate their 
susceptibilities \. 



t=i 



(5) 



here and in the following, V is the spatial lattice volume. 
The chiral condensate vp-ip is given by[l2j]: 



I 4nN f I 1 
-Re ( Tr 

VN t 2 \ Mt 



(6) 



The susceptibility of Polyakov loop xl is: 

XL = V((L-(L)f), (7) 

The susceptibility of plaquette variable \p and the sus- 
ceptibility of chiral condensate x^ip are defined as : 



XP = VN t ((P~(P)) 2 ), 
= VN t {{W-{W)) 2 ) 



(8) 
(9) 



where Nt is the number of temporal sites of the lattice. 

At high temperature, QCD with massless quarks is 
believed to restore the chiral symmetry which is spon- 
taneously broken. This is the chiral transition and the 
chiral condensate is the order parameter. However, due 
to the fact that our definition of chiral condensate for 
Wilson fermions is the naive definition and the Wilson 
fermions explicitly breaks the chiral symmetry, the mean- 
ing of (ipip) at k ^ k c is not clear. One should make a 
subtraction to compensate for the additive renormaliza- 
tion of the quark mass[l3l|. A properly subtracted chi- 
ral condensate (ipip) can he defined via an axial vector 
Ward-Takahashi identity[14]. Nevertheless, we employ 
the naively defined (?pip) an d the susceptibility on which 
we don't make a subtraction to compensate for the influ- 
ence of the Wilson term. 

However, when the system is at crossover or criticality, 
these physical observables will display sharp changes and 
their susceptibilities will display a peak, from which we 
determine the transition point. 

In a finite volume, the susceptibilities are always an- 
alytic functions, even in the regime where phase transi- 
tions occur, however, in the infinite volume limit, phase 
transitions reveal themselves through the divergences of 
the susceptibilities, whereas for crossover, susceptibili- 
ties are finite. The order of the transitions can be de- 
termined by the finite size scaling of the susceptibilities. 
The susceptibility at transition point Xmax behaves as 
Xmax oc V a , with a the critical exponent. If a = 0, the 
transition is just a crossover; If < a < 1, it is a second 
order phase transition; If a = 1, it is a first order phase 
transition, accompanied by the double peak structure in 
the histogram of the observable and flip-flops between 
the two states in the MC history [l5l]. 

Since the effect of the Wilson term is not subtracted 
and its volume dependence is non-trivial, whether the 
finite volume scaling behavior of the chiral susceptibility 
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is consistent with the scaling behavior described above is 
an open question. 

We also calculate the chiral condensate which is defined 
via an axial Ward-Takahashi identity [3], we will refer 
to it as subtracted chiral condensate and denote it by 
(4>il>) sub, and this properly defined (ipip)sub was employed 
in Ref. pi, 03, 



(10) 



where Z is the normalization coefficient, and the tree 
value of it is (2k) 2 which is sufficient for our study. The 

"ITT 



(ii) 



current quark mass m q is defined through [F 

2m q (0\P\n(p = 0)) - -m T (0|A 4 K(p - 0)>, 



where P is the pseudoscalar density tpj^ip, A4 is the 
fourth component of the local axial vector current 
■07574'i/), and \tt) and |0) stand for the pion and vacuum 
state, respectively. On the lattice, 



with 



R(z) 



2m q — lim R{z), 



(E»,„,t^(s,y,z,*MQ)) 



(12) 



(13) 



III. MC SIMULATION RESULTS 

In this section, we will present our results for simulat- 
ing QCD with two degenerate flavors of Wilson fermions 
at finite temperature T and imaginary chemical poten- 
tial ifxi. The HMC algorithm is used(20|. To deter- 
mine the pseudo-transition point /3 c (a/ijV, we use the 
Ferrenberg-Swendsen reweighting method|19|. The sim- 
ulations are performed on the V x N t = 8 3 x 4 lattice 
at K = 0, 0.005, 0.165. The molecular dynamics time 
step St is chosen in such a way that the acceptance 
rate is approximately 80 — 90% otherwise stated. There 
are 20 molecular steps for each trajectory. We gener- 
ate 20,000 trajectories after 5,000 trajectories of warmup. 
Ten or twenty trajectories are carried out between mea- 
surements. To determine the order of phase transition 
at some parameters, larger lattices are also used for fi- 
nite size scaling. When calculating the quark mass m q , 
we perform simulations on the 8 2 x 20 x 4 lattice while 
keeping other parameters unchanged. We use the con- 
jugate gradient method to evaluate the fermion matrix 
inversion. 



A. RW TRANSITION AT IMAGINARY 
CHEMICAL POTENTIAL 

In this section, we present the results of simulation 
for addressing the Z(3) transition, and the simulation is 
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FIG. 1: Histogram of tp/n at RW transition point a\ii 
7r/12 ps 0.262, where if is the Polyakov loop phase. 
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FIG. 2: (</?) /n as a function of a\xi for some different values 
of 8. 



performed with St = 0.02 for which the acceptance rate 
is approximately 90 — 95%. 

The SU(3) gauge theory with fermions at imaginary \x 
has periodicity with period 27rT/jV c ||, 0,111]. In the m S n 
temperature deconfined phase, there is a first order phase 
transition between different Z(3) sectors, while in the low 
temperature phase, the transition becomes a crossover 
at some critical imaginary chemical potential values fi c j 



EL 

T 



2tt 



k 



(14) 



The different Z(3) sectors can be distinguished from each 
other by the phase of Polyakov loop. In our case, i.e., 
N c = 3 and N t = 4, the first Roberge- Weiss (RW) tran- 
sition to different Z(3) sectors should appear at a/// = 
7r/12 w 0.262. Because the system will tunnel into the 
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unphysical Z(3) sector above pi/T = 7r/3, our method is 
limited up to pj/T = tt/3. 

Fig. [1] shows the history and probability distribution 
of the phase ip of the Polyakov loop at api = 0.262, k = 
0.165 and (3 — 5.260. Figure shows i^/i as a function 
of api at some different values of j3. These indicate that 
at afij w 0.262, and T > Te (where f3 is larger than 
[5.245, 5.255] ), there is a first order phase transition with 
(ip) changing between values of and — 27r/3. 



B. DECONFINEMENT TRANSITION AT 
IMAGINARY CHEMICAL POTENTIAL 

In order to investigate the deconfinement transition, 
we take measurements of plaquette variable P, Polyakov 
loop norm \L\, chiral condensate i/tip and their suscepti- 
bilities xp> X\l\, Xm-> an( i the subtracted chiral conden- 
sate {tpip) su b in the first Z(3) sector api < (n/3N t ) at 
api = 0.0, 0.10, 0.14, 0.18, 0.21, 0.24 at k = 0.165 by 
using the Ferrenberg-Swendsen reweighting method. 

The values of (3 at which we make simulations for the 
Ferrenberg-Swendsen reweighting method and the quark 
and meson 7T, p screening mass are presented in Table. U 
except for api = 0.10, 0.18. In order to calculate the 
subtracted chiral condensate, we must know the quark 
mass first. At api = 0.10, 0.18, the values of (3 at which 
we make simulations are the same as those at api — 
0.0 and we use the quark masses obtained at the four 
different values of (3 at a pi — 0.0 and a pi = 0.21 as the 
quark masses at api = 0.10 and api — 0.18, respectively. 
For the quark mass differs slightly at the same (3 and 
different api. 

The values of plaquette, Polyakov loop norm, chiral 
condensate and their susceptibilities are plotted in Fig. [3] 
and Fig. [H respectively (we only plot them for three val- 
ues of api for clarity). We also display the values of the 
subtracted chiral condensate for only three values of api 
in Fig. [5] These observables at other api's have similar 
behavior as shown in Fig. 1314151 

From Fig. [3] and Fig. O one sees that around the same 
/3's, the values of P, \L\, •ipip and the subtracted chiral 
condensate {ipip) a ub change rapidly and the value of ipip 
is larger than that of (ipi()) su b at the same (3 and api. 
Fig. |4] tells that the locations of the peaks for %Pi X\lU 
Xm are consistent with each other within errors. We de- 
termine the transition points (3 c {api) from the locations 
of susceptibility peaks, the results are listed in Table. HT1 

In Ref. @, it has been established in detail that be- 
cause the partition function Z as an even function of 
ap leads to an even susceptibility x, and at the transi- 
tion points, dx/9/3 = 0, this expression implicitly defines 
(3 c {ap) as an even function of the real chemical poten- 
tial ap due to implicit function theorem; that when the 
purely imaginary chemical potential is considered, the 
considerations are unchanged, so, pseudo-critical line of 
the transition at imaginary chemical potential is simply 
the analytic continuation of the pseudo-critical line at 



TABLE I: Results of 71", p meson and twice quark screening 
mass for Nf — 2 on 8 2 x 20 x 4 lattice. The acceptance rates 
are approximately 75 — 83%, with the exception at a/ii — 
0.24, /3 — 5.24, at that point, the acceptance rate is 70%. 



api 


P 




2m q a 


m p a 


0.00 


5.195 


1.244(2) 


0.4193(9) 


1.354(2) 




5.215 


1.301(2) 


0.2550(5) 


1.410(1) 




5.235 


1.327(2) 


0.1893(4) 


1.437(1) 




5.255 


1.354(2) 


0.1386(2) 


1.455(1) 


0.14 


5.195 


1.203(2) 


0.4542(9) 


1.298(2) 




5.215 


1.217(2) 


0.3286(7) 


1.290(2) 




5.235 


1.301(2) 


0.2024(3) 


1.328(2) 




5.255 


1.320(2) 


1546C2) 


1 345(2) 


0.21 


5.195 


1.182(3) 


0.461(1) 


1.274(2) 




5.215 


1.140(2) 


0.421(1) 


1.228(2) 




5.235 


1.177(3) 


0.2527(7) 


1.201(2) 




5.255 


1.278(2) 


0.1558(3) 


1.237(2) 


0.24 


5.200 


1.169(2) 


0.456(1) 


1.263(2) 




5.220 


1.117(3) 


0.416(1) 


1.211(2) 




5.240 


1.098(3) 


0.338(1) 


1.148(3) 




5.260 


1.242(3) 


0.1013(2) 


1.192(2) 



real chemical potential; hence that the pseudo-critical 
transition line at imaginary chemical potential (3 c (api) 
is an even function of api and can be fitted well by a 
polynomial of degree one in (api) 2 without taking into 
account the term of degree two in (api) 2 , that is to say: 

Pc{api) = c + ci(api) 2 + 0(a 4 p 4 I ) 1 (15) 

After we obtain the expression for (3 c {api) as a polyno- 
mial of [api) 2 , we continue back to the real chemical 
potential and get /3 c (ap) as a function of ap. 

We use the least squares method to fit the data in Table 
im the coefficients and \ 2 /dof are listed in Table [iTTl The 
fitting range and the line are presented in Fig. [Sj From 
Fig. O we find that the coefficients of terms in (api) 2 
with higher order than one are difficult to be determined 
with high precision. From Table IIII1 we can find that 
the fitting result from chiral condensate is better than 
the results from P and L, so our choice for the pseudo- 
critical transition line is: 

[3 c {api) = 5.201(3) + 0.722(80)(a^/) 2 + 0{a 4 pj). (16) 

the errors are the fit errors. 

We estimate the pseudo-scalar meson mass m w , the 
vector meson mass m p and their ratio m w /m p at our 
simulation points from the data in Ref. 21]. By using 
the standard quark and gauge action, Bitar et al. found 
that at k — 0.16, (3 = 5.28, om, = 1.213 ± 0.004 and 
am p = 1.287 ± 0.0005, at k = 0.17, (3 = 5.12, am^ = 
1.088 ± 0.003 and am p = 1.210 ± 0.005. Our critical 
(3 values range from (3 = 5.199 to [3 — 5.243 at n = 
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FIG. 4: Susceptibilities for the plaquette, the Polyakov loop 
FIG. 3: Mean values of the plaquette, the Polyakov loop norm norm and the chiral condensate at k = 0.165. 
and the chiral condensate at n — 0.165. 



0.165, so we estimate that at the transition points in our 
simulation, am,'s are in the interval from 1.213 ± 0.004 
to 1.088 ± 0.003, amp's from 1.287 ± 0.0005 to 1.210 ± 
0.005, and the ratios of m^/nip are between 0.943(3) and 
0.899(4). 



C. DECONFINEMENT TRANSITION AT REAL 
CHEMICAL POTENTIAL n 

Now it is trivial to get the pseudo-critical line on the 
(p,T) plane. Because /3 c (a/ij) is an analytic function of 
a/z/[6j, we can analytically continue from the imaginary 
chemical potential to the real one. Replacing [i\ by — i/j, 
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TABLE II: Collection of pseudo-critical transition points at k = 0.165, determined by locating the peaks of the susceptibilities. 





0.00 


0.10 


0.14 


0.18 


0.21 


0.24 


Pc from xp 


5.199(5) 


5.208(4) 


5.215(7) 


5.226(4) 


5.233(3) 


5.243(3) 


Pc from xl 


5.199(5) 


5.209(4) 


5.215(5) 


5.226(4) 


5.233(3) 


5.243(3) 


Pc from Xi>4, 


5.200(5) 


5.208(4) 


5.218(7) 


5.226(5) 


5.234(4) 


5.242(3) 



To translate our result into physical unit, we use the two 
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FIG. 5: Mean values of the subtracted chiral condensate at 
k = 0.165. 

FIG. 7: Illustrative figure of transition temperature as a func- 
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m Eq. we obtam /3 c (a M ), = (19) 

Pc(a/j) = c - ci(a[i) 2 + 0(a 4 fi 4 ) 

= 5.201(3) - 0.722(80)(a/i) 2 + 0(a 4 /x 4 ).(17) by replacing a with l/(N t T c ), it gives: 
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Tc(0) 



exp 



1 61 . / 1 



exp 



126 r 



2bl 



In 



1 

66^ 



(20) 



Expand the right hand side of Eq. ([20]) as a series of fi 2 
and neglect the higher order terms , we can obtain the 
expression of T c as a function of ji 2 B . 

T T ( cUiB \s =l- 0.00722(80)(^) 2 , (21) 
TcIhb = 0) T 

where the baryon chemical potential [Ib is related to the 
quark chemical potential by jj, = hb/N c and the error 
only reflects the error on c\. T c (m3 — 0) is set by the 
critical temperature for 2-flavor QCD at ms = 0. 



TABLE III: Coefficients of Eq. flgjl by fitting the data in 
Table. HU 





CO 


Cl 


X 2 /dof 


/3c from \P 


5.200(3) 


0.748(79) 


0.057 


Pc from xl 


5.201(3) 


0.738(78) 


0.075 


(3 C from x^ 


5.201(3) 


0.722(80) 


0.104 



Recently, Bernard et al. [23J studied the transition tem- 
perature of 3-flavor, (2+1) flavor QCD, they obtained 
T c = 169(12)(4)MeV or T c = 174(ll)(4)MeV for (2+1) 
flavor. Cheng et al.[24| performed calculation of the 
transition temperature of (2+1) flavor QCD, and T c — 
192(7) (4)MeV is their result. Karsch, Laermann and 
Peikert obtained T c (^ = 0) = 173(8)MeV in the chiral 
limit for staggered fermions[25|. Ali Khan et al. used a 
renormalization group improved gauge action and clover- 
improved Wilson quark action to investigate 2-flavor 
QCD and they obtained TJji = 0) = 171(4)MeV[]J]. 
The result of Karsch et al. and Ali Khan et al. are 
consistent with each other. If we take 173(8)MeV as 
T c {lJ<B = 0), then the transition temperature for Nf = 2 
is described by a line illustratively plotted in Fig. [7] from 
which we can see that T c decreases with increasing fig. 

The imaginary chemical potential method is valid in 
the range hb/T < 7T, and the pseudo-critical (3 is a poly- 
nomial of (a/i/) 2 . However, the data in Fig. [6] imply 
that we can only calculate the first two coefficients of the 
polynomial with high precision, so the continuation from 
imaginary chemical potential to real one is restricted in 
the range of small an and therefore Eq. (f2Tj) is valid in 
the range of small fiB- 



D. NATURE OF PHASE TRANSITION 

In order to determine the nature of the phase transition 
with imaginary chemical potential, we investigate the his- 
tory, histogram and finite size scaling of MC simulation 
at K = 0, 0.005, 0.165. On lattice 8 3 x 4, 12 3 x 4, 16 3 x 4, 
at k = which corresponds to quenched limit or pure 
gauge theory, we find the critical /3 is 5.70 by determin- 
ing the location of the peak of X\ LL The value of /3 C is 
consistent with the result of Ref. |26l . I27T ] . We plot the 
history and histogram of (\L\) around critical (3 C — 5.70 
in Fig. [51 From the histogram and MC history of {\L\), 
we see that near /3 = 5.70, there are two-state signals 
which arc an indication of first order phase transition. 

Because at k = 0.0, quarks have no effect on the 
system, it is natural that the value of (3 C and the two- 
state signal are the same for other values of afii in the 
quenched limit. So we conclude that at other values of 
a/i/ in the quenched limit, the phase transition is of first 
order. We also make simulations at k = 0.005 on lattice 
8 3 x 4, the result is presented in Fig. [9] which tells us that 
for very heavy quarks, the system with imaginary chem- 
ical potential has the feature of first order transition. 

In order to estimate the lattice spacing at (3 — 5.70, 
we use the results in Ref. [27j from which we know that 
at /3 = 5.6925(2) with the temporal extent N t — 4 at the 
infinite volume, y/aa = 0.4179(24), a is the string ten- 
sion. Using-yo- = 420MeV, we estimate that the lattice 
spacing aw 0.995G6V- 1 . 

At k = 0.165, am = 0.24, the critical is 5.242(3) 
or 5.243(3), and (3 — 5.244 is consistent with them 
within errors, so at ^ = 5.244, we can evaluate 
the spatial dependence of susceptibilities of Polyakov 
loop norm and its time history and histogram on 
lattice of spatial size of 8 3 , 10 3 , 12 3 , 14 3 , 16 3 , 18 3 , 20 3 
with temporal extent Nt — 4 with acceptance rate 
82%, 81%, 75%, 68%, 67%, 55%, 51%. We generate 
1,000 configurations except on lattice 12 3 x 4 where 2,000 
configurations are generated. We present the spatial de- 
pendence in Fig. |TD] and the time history and histogram 
on lattice 8 3 x 4, 12 3 x 4, 16 3 x 4 and 20 3 x 4 in Fig. QTJ 

From Fig. [TUlwe find that the peak heights of Polyakov 
loop norm susceptibilities are approximately the same ex- 
cept for spatial volume 12 3 , 18 3 . On lattice 12 3 x 4, as a 
comparison, we measure after the first 1,000 config- 
urations are produced, we find that — 21.74(9.46), 
while x\l\ = 16.38(5.47) after the statistic is doubled. We 
can expect that when statistic is large enough, the val- 
ues of x\l\ an d their errors on lattice 12 3 x 4, 18 3 x 4 will 
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FIG. 10: Peak heights of susceptibility of Polyakov loop norm 
as a function of spatial volume at a\ii = 0.24. 



FIG. 8: Time history and histogram of Polyakov loop at k 
0.0. 
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FIG. 9: Time history and histogram of Polyakov loop norm 
at k = 0.005. 



decrease. The history and histogram of Polyakov loop 
norm plotted in Fig. 111! together with the peak height 
change with spatial volume, shows that the transition at 
a/i/ = 0.24 is a crossover. 



IV. DISCUSSION 
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We have studied the phase diagram of lattice QCD 
with the two flavor of Wilson fermions through the sim- 
ulations with imaginary chemical potential. In this case 
the partition function is periodic in imaginary chemical 
potential. The different Z(3) sectors are characterized by 
the phase of Polyakov loop. The Z(3) transition which 



FIG. 11: History and histogram of Polyakov loop norm at 
different spatial volume at afii = 0.24. 
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occurs at fXi/T = 2ir(k + l/2)/3 is of first order in the 
high temperature phase. 

Our study shows that there is a first order phase transi- 
tion at k = which corresponds to infinite heavy quarks 
or the quenched limit, it is natural that the (3 C and 
hence the critical temperature have no dependence on 
the chemical potential based on the fact that the fermions 
have no effect on the system in the quenched limit. 

From the experience and literature, we expect that in 
general, the lighter the quark mass, the stronger effect 
on the system the fermions have. At At = 0.165, we ob- 
serve that the values of P, \L\, 't/jip and subtracted chiral 
condensate (tpip} S ub change rapidly around f3 c and the 
transition points determined from the susceptibilities of 
P, \L\, ipip coincide with each other within errors. The 
transition at k = 0.165 which corresponds to a value of 
ratio of m^/nip in the range from 0.943(3) to 0.899(4) at 
imaginary chemical potential a\ii — 0.240 is a crossover, 
as discussed in the preceding section. 

As for the transition temperature as a function of 
chemical potential, as discussed in Ref. c 6j, the critical 
line can be well described by a linear function of [i 2 B . We 
make simulation at k = 0.165 and have not investigated 
the quark mass dependence of our results. Our central 
result is represented by Eq. (|2"Tj) which is qualitatively 
consistent with yet quantitatively slightly different from 
that in Ref. Q taking errors into account. We think 
that it is probably because our simulation is at a point of 
quark mass larger than that in Ref. @. At our simula- 
tion points, the ratio of pseudo-scalar meson mass to 
vector meson mass m p is between 0.943(3) and 0.899(4), 
these large ratios mean that the quark mass is large at 
our simulation points. 

In the process of deriving Eq. (|2"Tj) , we use the 2-loop 
perturbative solution to the renormalization group equa- 
tion between lattice spacing a and (3. However, in our 
simulations on lattice with N t = 4, the values of critical (3 
range from 5.199 to 5.243, the coupling is so strong that 
using the 2-loop expression is not a good choice. One 
would need the non-perturbative expression between o 
and P, but it is not determined so far. So we have no 
choice but to use the 2-loop perturbative expression be- 

1 One of the authors Wu thanks Philippe de Forcrand for telling 
Wu those information. 



tween a and j3. It is known qualitatively that the lattice 
spacing varies faster than predicted by the 2-loop per- 
turbative formula at strong couplings. This will have the 
effect of increasing the curvature in Eq. ({2T]) 1 . 

Solving Eq. (f2"Tj) , we can obtain the T(^b) as a func- 
tion of [i B . We take 173(8)MeV as T C (^ B = 0) and illus- 
tratively plot the transition temperature T c (MeV) ver- 
sus baryon chemical potential /^(MeV) in Fig. [7] from 
which we find that the transition temperature T c de- 
creases slowly with increasing (is- This behavior is in 
accordance with the physical picture. With baryon den- 
sity increasing, the interaction between quarks and glu- 
ons becomes weaker and thus quark and gluon degrees 
of freedom get more easily excited, therefore, the critical 
temperature decreases with increasing baryon chemical 
potential. As discussed in Sec. IIII C[ Eq. (|21[) is valid in 
the range of small /j,b- 

In order to get the transition occur for small m^/nip 
and make the use of 2-loop perturbative relation between 
lattice spacing a and (3 more reliable, lattices with larger 
temporal extent iV t would be used. The investigation of 
chemical potential dependence of transition temperature 
in the chiral limit and quark mass dependence of tran- 
sition temperature awaits further work. Moreover, how 
to extract the information about the nature of transi- 
tion with real chemical potential from the behavior with 
imaginary chemical potential remains an open question. 
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